Universal Velocity Profile for Coherent Vortices in Two-Dimensional Turbulence 
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Two-dimensional turbulence generated in a finite box produces large-scale coherent vortices co- 
existing with small-scale fluctuations. We present a rigorous theory explaining the = 1/4 scaling 
in the V oc r~^ law of the velocity spatial profile within a vortex, where r is the distance from 
the vortex center. This scaling, consistent with earlier numerical and laboratory measurements, 
is universal in its independence of details of the small-scale injection of turbulent fluctuations and 
details of the shape of the box. 
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The generation of large-scale motions from small-scale 
ones is a remarkable property of 2D turbulence. This 
phenomenon is a consequence of the energy transfer to 
large scales [H, 0, realized via inverse cascade. Simula- 
tions 0, |3l and experiments 0, 01 show that accumulation 
of energy in a large-scale coherent structure is observed 
at sufficiently long times if the friction is small enough 
and docs not prohibit the energy cascade from reaching 
the system size. Recent interest in undcrstandingthe 
structure of this state was sparked by experimental 
and numerical observations of large-scale coherent 
vortices associated with energy condensation in forced, 
bounded flows. One motivation for studying 2D tur- 
bulence comes from its structural and phenomenological 
similarity to quasi-geostrophic turbulence lllUl2 |. such 
as that observed in planetary atmospheres [l3|. Also as 
suggested in 0, the emergence of large-scale coherent 
structures in 2D is related to the confinement transition 
in magnetic plasmas whose slow dynamics is described 
by quasi-2D hydrodynamic equations. 

In this Letter, we examine the large-scale vortices, gen- 
erated by inverse energy cascade in a finite box. We 
begin our discussion with a brief review of the classi- 
cal theory of inverse cascade by Kraichnan, Leith, and 
Batchelor (KLB) P,[3,[il. The essential difference of 2D 
turbulence and 3D turbulence is the presence in the for- 
mer of a second inviscid quadratic invariant, in addition 
to energy, the enstrophy. Therefore, stirring of 2D flow 
generates an enstrophy cascade from the forcing scale, I, 
down to smaller scales (direct cascade) and also gener- 
ates an energy cascade from the forcing scale up scales 
(inverse cascade). Viscosity dissipates enstrophy at the 
Kolmogorov scale, rvisc, which is much smaller than I 
when the Reynolds number is large. In an infinite sys- 
tem, the energy cascade is eventually blocked at the scale 
of Tfric by friction, thus resulted in establishing the two- 
cascade stationary KLB turbulence for rfric ^ I- The 
Kolmogorov phenomenology (see, e.g., for a review), 
KLB predicts the velocity power spectrum in the di- 
rect cascade and k'^^^ in the inverse cascade, where k is 
the modulus of the wave vector. These KLB theoretical 



predictions were confirmed in simulations |15j and labo- 
ratory experiments [3, 17, 3]. (Note also discussion of 
experimental evidence of simultaneous inverse and for- 
ward cascade in the infinite system setting [i3|-) If the 
frictional dissipation is weak and rfric exceeds the system 
size L, then ultimately the KLB regime is not applica- 
ble and a large-scale coherent flow (occasionally called a 
condensate) emerges (ioj . 
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FIG. 1: Average vorticity profile observed in simulations (A, 
Fig. 2C from [l3|) and experiment (B, Fig. 3 from [^). 

Laboratory experiments showed that the coherent flow 
contains one or two vortices, depending on the boundary 
conditions, and taking roughly a half of the system size [9|, 
21 1 . Numerical simulations of [l3| reported a well-defined 
scaling for the average velocity profile in the interior, V oc 
r^^ , as a function of distance, r, from the vortex center, 
and rj w 1/4. Similar scaling was also observed in the 
thin layer experiment Fig. [1] summarizes the results 
of [l,[l3| for the average vorticity, CI, demonstrating the oc 
7.-5/4 bgjiavior corresponding to rj ~ 1/4. (We mention 
the experimental results of [9| to emphasize emergence 
of the scaling range, possibly not reached in preceding 
experiments, e.g. |2l| . presumably because of a somehow 
higher surface friction. ) 

Our main result is a rigorous derivation of the 77 = 1/4 
scaling and explanation of why this scaling is universal. 
Our derivation is the first of the kind in the field of tur- 
bulence as predicting universal scaling for a structure 
emerging as the result of a nonlinear balance between the 
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small-scale turbulence and the coherent structure gener- 
ated by turbulence. The key feature that allowed us to 
derive this analytical result, is the smallness of the am- 
plitude of the background velocity fluctuations in com- 
parison with the coherent part. In essence, this small 
parameter provides an asymptotically accurate trunca- 
tion of the generally infinite system of Hopf equations 
on the level of third order correlation functions of veloc- 
ity /vorticity. The 77 = 1/4 scaling emerges from an ex- 
plicit solution of the resulting system of equations. The 
main contribution to the third-order Hopf equation (for 
the third order correlation function) is associated with a 
zero mode of the respective integro-differential operator 
representing the homogeneous part of the equation. This 
result is substituted into the second-order Hopf equation, 
thus treated as a linear inhomogeneous equation for the 
pair-correlation function, with the third-order correlation 
function calculated on the previous step. Similarly, the 
first order Hopf equation is a linear inhomogeneous equa- 
tion with respect to the first moment (the coherent term), 
resulting in a closed expression for the scaling exponent. 
Our strategy below is to derive the set of equations, in- 
troduce the truncations and show that a scale-invariant 
expression with rj = 1/4 gives a solution consistent with 
the assumptions made in the process of the derivations. 

The 2D velocity field v is assumed controlled by the 
Navier-Stocks equation 



9fti> + {vV)lj — —auj + vV'^uj + curl /, 



(1) 



formulated in terms of the vorticity lo — cwclv. One as- 
sumes that the fluid is incompressible, S/ ■ v = 0. The 
terms on the right-hand side of Eq. ([T]) represent the 
bottom friction, the viscosity and the turbulent forcing 
respectively. The force per unit mass, /, is assumed 
to be random, zero mean, statistically homogeneous in 
space and time, and correlated at an intermediate scale, 
I, called pumping scale. We study the case in which the 
pumping scale is much smaller than the size of the sys- 
tem (of the box), L, and it is much larger than the Kol- 
mogorov (viscous) scale, rvisc- The energy density (per 
unit mass) e injected by the forcing / in a unit of time 
per unit mass is considered constant in space and time. 

If one starts from zero velocity and turns on the pump- 
ing / at t = 0, in time ^ P^^/e^^^ a direct enstrophy 
cascade is established in the range of scales between the 
pumping scale and the viscous scale. The establishment 
of this direct cascade is followed by a much slower growth 
of the inverse energy cascade from the pumping scale to 
larger scales. The energy-containing scale of the inverse 
cascade grows, as knv ~ e^^^t^^^, until the scale reaches 
the system size L at time Tmh ^ {L^/eY^'^. After that 
the system has no choice but to deviate from spatially 
homogeneous KLB regime, producing a large-scale co- 
herent flow. This picture is correct provided the bottom 
friction is sufficiently weak, a <^ e^^^L~^/^, as assumed 
in the following. Establishing the spatial profile of the 



resulting average velocity at times t > Tinu is our main 
objective. 

Once the large-scale coherent flow has emerged, the 
total velocity field, v{t,r), can be decomposed into co- 
herent and fiuctuating parts, v{t,r) = V{t,r) + u{t,r). 
By definition of the coherent part, (u) = 0, where an- 
gular brackets indicate averaging over the temporal scale 
of the vortex turnover time, ~ L/V. After a transient 
period, i.e. once the large-scale flow has matured, the 
injected energy is mainly accumulated in the coherent 
component of the velocity V{t, r), which grows slowly in 
time, oc -s/i, corresponding to the linear in time growth of 
the total energy. Ultimately, the average velocity profile 
V is stabilized by the bottom friction, then the velocity 
amplitude is determined by the balance between energy 
injection and dissipation, where thus V oc Xj ^/a. Let us 
consider averaging Eq. ([T]) 



(2) 



where 17 = (cj) is the average vorticity and vj = curl it 
stands for respective fiuctuations. In Eq. ^ we have 
ignored viscosity and dropped pumping, both irrelevant 
at scales larger than I. Our description, detailed below, 
is based on the assumption that the coherent flow dom- 
inates fluctuations, V u, possibly everywhere except 
for a small neighborhood of the vortex core. (This as- 
sumption can be accurately verified via a self-consistent 
multi-step procedure including an analysis of the multi- 
point correlation functions. A detailed discussion of these 
technical details is postponed for a longer publication.) 

In the periodic set-up, e.g. realized in simulations [lO|, 
a pair of vortexes forming a dipole is formed, whereas 
in a bounded box set-up one observes a single vortex, 
as seen in the experiment of 0]. Other structures, e.g. 
more than one pair of vortices or stripes or combina- 
tions of stripes and vortexes, can also emerge in boxes 
of special shapes, such as those with large aspect ratios 



or non-trivial topology (e.g., stripes and rings) |22l . 
In the following we will focus on an analysis of a vortex 
which applies equally well to either of the two cases men- 
tioned above. (Note also that our approach to describing 
the shape of the vortex is based on analysis of stochastic 
Navier-Stokes, and as such it is distinctly different from 
the quasi-equilibrium 2d approach 2^ 2J, , postulat- 
ing a distribution of Gibbs-kind controlled by the set of 
lagrangian multipliers associated with different moments 
of vorticity.) 

An emergence of the coherent vortices results in an in- 
homogeneous redistribution of energy. After a vortex (or 
pair of vortexes) has emerged, the global profile of V{t, r) 
shows two regions, corresponding to a vortex exterior and 
an interior respectively. In the vortex exterior the average 
velocity is estimated as y ^ Vet orV^ \f^Jo<- (the latter 
corresponds to the stationary case, in which the turbu- 
lence is stabilized by friction), while inside the vortex 



3 



the coherent part is much larger and (up to small varia- 
tions wc are ignoring) its components are = V(r) and 
Vr = 0. Inspired by the results of the numerical and 
laboratory [9[ experiments, we assume that the spatial 
profile of the coherent part in the interior of the vortex 
is algebraic, that is 



V{t,r) = Vo{L/r)\ 



(3) 



where the distance r is measured from the vortex center 
and Vo estimates the coherent part of the velocity in the 
vortex exterior. Eq. ([3]) is correct for the rcore ^ r <C i 
range. Here Tcoio is the size of the vortex core. The 
term, in Eq. Q is zero due to the isotropy of the 

vortex. Therefore the vortex profile is determined by a 
balance of the first and third terms in Eq. ([2|) . Obviously, 
Eq. ^ is not closed and one naturally needs to consider 
an additional equation for the pair correlation function 
of velocity /vorticity fiuctuations inside the vortex. 

In fact, it is convenient to derive these extra equations 
for the averages in two steps, first rewriting Eq. ([1]) 

{dt + a)ur + d^N-^tur - d^r^M-\uVuj) = , (4) 

where both the force and the viscosity terms are dropped. 
The differential operators in Eq. ^ are 

M^r[id, + ir+dl], (5) 

ic = v{dl + 2d, + 2 + dl)-{dlv), (6) 

where g = ln(r/L). Then, wc introduce the pair correla- 
tion function of the radial velocity fluctuations 

<^{t, ri, r2, </?) = {ur{t, ri, Lpi)ur{t, ra, </J2)) , (7) 

where = Lp\ — ip^. The correlation function is invariant 
under the transformation, — > — i^, r\ *^ r2, correspond- 
ing to the permutation of the points labeled by 1 and 2. 
Using this property and assuming analyticity of the pair 
correlation function ([7]) for small ip and p = ln(ri/r2), we 
derive the following expression for the single-point cross- 
object of the second-order, (urvo), appearing in Eq. ([2|): 



(8) 



^=0,(^=0 



where r = ^T\r2 and J5 = 1 -|- rdrjl. Note that only 
antisymmetric in p> term in $ contributes to {urw). 

Multiplying Eq. ([4]) by the velocity at another spatial 
point and averaging the resulting equation over fluctua- 
tions, one derives 

U^^U^^ (aAs-^i -M/C2) $(ri,r2,(^) 

= rjTV'f ^Vi {u{rx)w{ri)ur{r2)) 
~rlN2^^2{ur{ri)u{r2)Tu{r2)) , (9) 

where the irrelevant (asymptotically small) terms, con- 
taining the time derivative and the friction coefficient a. 



are omitted. The operator on the left hand side of Eq. ([9]) 
can be rewritten as 

N2ti - M1IC2 = {r2Vi - riV2) X 



-AB^dl + 25(1 - r]^) coth[(l + r])p/2]dpj. (10) 

When the separation ri — r2 is sufficiently small, the 
right hand side of Eq. ([9]) controls the inverse energy 
fiux, exactly as in the traditional KLB case. Indeed, 
in the spatially homogeneous case, the correlation func- 
tion {ua{ri)Tu{ri)up{r2)) ~ e depends solely on ri — r2 
and it is also divergentless due to Vit = 0. Substitut- 
ing the Kolmogorov estimate $ — ^(O) ~ — '"2!)^^'^ 
into the left hand side of Eq. ([9]), one finds that the 
term is negligible for \ri — r2\ = ri2 C*i where 
(i,{r) ~ Q^l^r'^l'^t'^l^iyT jVy"^!"^ is thus an important 
scale dependent on the distance from the vortex core, r. 
One concludes that for ryi ^ ^-k(y)^ the inverse cascade 
is modified by the coherent fiow. 

However, due to isotropy the term in $ related to the 
KLB cascade does not contribute to the object of our 
prime attention, {uvj). First, wc look for such solutions, 
also remaining regular at small ri2, in terms of the zero 
modes of the operator on the left-hand side of Eq. ([5]), 
thus ignoring the smaller right-hand side in the equation. 
However, these zero modes do not contribute to (urw), 
because the last term in the operator PU)) prohibits odd 
in p zero modes to be regular for small p. Therefore, 
to extract a non-zero contribution to {urw), one has to 
account for a correction to $, (5<I>, related to the right- 
hand side of Eq. ([9]). To get a non-trivial contribution 
one ought to carry our analysis to the next order in the 
Hopf hierarchy describing the triple velocity correlation 
function, F (fr(ri)w,.(r2)u,.(r3)). 

The principal terms in the third-order Hopf equation 
governing F are 



-^N^'t, + ^N2't2 + ^M^'tA F 
dpi dp2 dp3 



0. 



(11) 

where we again omitted asymptotically irrelevant terms, 
including the time derivative term, the friction term and 
also the contribution related to the fourth order correla- 
tion function. Formally, any zero mode of the operator 
K satisfies Eq. ([TT|) . and the quest is to find the scale- 
invariant zero mode of K, = exp^irmp + (3mQ) where 
m is an integer (to guarantee smoothness at the smallest 
scales) and f3m ~ yjvn? +Tf — 1 — 1, which generates a 
nonzero contribution into {urVj) via Eqs. (j9l8p and has 
the smallest possible . The first terms in the hierarchy 
of possible candidates are 



F (X Zmiri, Pl)Zk{r2, P2)Z-rn-k{r3, ips) + 



(12) 
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where the dots represent the sum of the terms that are 
obtained from the first product in Eq. (|12p by permuting 
the indices 1,2,3. However, the expression generates 
an odd outcome for the right-hand side of Eq. , thus 
resulting in an even correction 6^ giving no contribution 
to {urVj). Therefore, one has to look for higher order 
terms in the hierarchy. One finds that the desired zero 
mode can can be constructed with the help of an auxil- 
iary object, Xm = exp[im(p + (/?,„ -I- 1 -I- rfjg], satisfying 
(A/'m)~^/Cm-'^m = ^m^m, wherc A,n are real numbers: 



\-a^k-mZm{ri,(pi)Zkir2,(p2)X-jn-kir3,V3) + ■ 



,(13) 



and the dots stand for the sum of terms which accounts 
for respective permutations. Eq. (fT3|) is a solution of 
Eq. (HIl), provided ammAm + UkkAk - a^rn^k{m + 
k)A-m~k ~ 0. Choosing m = fc = 1 we find a term giv- 
ing a non-odd oc /•3')+v3+^-3 contribution to the right 
hand side of Eq. ([9]) . Then the correction to the pair cor- 
relation function is (5$ oc r^n+\^3+rf-2 ^ This result, fi- 
nally substituted into the last term in Eq. ([2]), translates 
into the 4?7 — 3 -I- + rj"^ = —rj relation, whose solution 
is ?7 = 1/4. To conclude, the ZZZ and ZZX terms, rep- 



resented by Eqs. (|T2|) and Eqs. (|T3|) . are the only struc- 
tures possibly contributing to the third-order correlation 
function F, and of the various allowed (nonzero) contri- 
butions, the ZZX term p3l) with m = k = 1 dominates 
F at ri,2 < L. 

Substituting the expressions, corresponding to i] ~ 
1/4, into the the Hopf equations of the first, second and 
third orders and estimating all the terms dropped in the 
derivation process confirms the validity of our asymptotic 
approximations. This completes our derivations. 

We now summarize our results. The main and some- 
what surprising result we just derived conserns universal- 
ity of the vortex mean profile. The scaling of the vortex 
shape is controlled primarily by a nontrivial zero mode 
of the operator on the right-hand side of Eq. pT|) and 
otherwise it follows from scaling relations between pairs 
of terms in the first- and second-order Hopf equations. 
Nothing in this solution is sensitive to the geometry of 
the box, or the details of the pumping. The solution also 
does not depend on the type of viscosity (hyper or nor- 
mal), or the damping coefRcient. Our conclusion does 
not depend of whether or not the coherent part grows in 
time or if it was already saturated by damping. Finally, 
our results make predictions going far beyond the main 
scaling statement, in particular detailed structure of an- 
gular harmonics is predicted for pair and triple correla- 
tion functions in the coherent regimes. Our theoretical 
statements call for accurate experimental and numerical 



tests. 

We conclude by mentioning a number of other compre- 
hensive questions raised by this study. Suppose a vortex 
or a pair of vortexes, internally tuned and built up from 
the energy flux are produced, but then the pumping is 
switched off. Will the initially formed vortex will keep 
its shape dynamically? Also if a somewhat different in 
shape, non-universal and large scale vortex is created, 
will it transform via decaying turbulence into the univer- 
sal shape predicted above? We conjecture that answers 
to both the questions are affirmative. These questions 
certainly require careful investigation in the future. 
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